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HIGHLIGHTS 


►  Mechanical  cell  testing  with  in-situ  short  circuit  location  identification. 

►  Macroscopic  jelly  roll  fracture  initiates  the  internal  short  circuits. 

►  Stress-based  fracture  and  short  circuit  criterion  developed. 

►  Finite  element  crash  simulation  of  the  cell  tests. 
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A  quasi-static  mechanical  abuse  test  program  on  cylindrical  Lithium  ion  battery  cells  has  been  performed 
at  a  state  of  charge  (SoC)  of  0%.  The  investigated  load  cases  involved  radial  crushing,  local  lateral 
indentation  and  global  three-point  bending  of  the  cell.  During  the  tests,  the  punch  load,  the  punch 
displacement,  the  cell  voltage  and  the  temperature  development  of  the  cell  have  been  monitored  using 
an  infrared  camera  and  temperature  sensors.  After  the  test,  the  cells  have  been  analysed  using  computer 
tomography. 

It  is  indicated  that  macroscopic  jelly  roll  fracture  on  a  global  scale  initiates  the  internal  short  circuits, 
revealed  by  a  sudden  decrease  of  the  global  mechanical  load  due  to  the  rupture,  followed  by  a  drop  of  the 
measured  voltage  and  immediate  increase  in  cell  temperature. 

A  macro-mechanical  finite  element  crash  simulation  model  has  been  established  for  the  cell  housing 
and  the  jelly  roll.  The  classical  stress-based  criterion  after  Mohr  and  Coulomb  (MC)  has  been  applied  to 
predict  fracture  and  the  initiation  of  an  internal  short  circuit  of  the  jelly  roll.  The  MC  criterion  correctly 
represents  the  punch  displacement  to  fracture,  where  the  predicted  fracture  locations  correspond  to  the 
observed  locations  of  the  internal  short  circuits  of  the  cells. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Plug-in  hybrid  electronic  vehicles  (PHEV),  range  extender 
battery  electronic  vehicles  (RE  BEV)  and  pure  battery  electric 
vehicles  (BEV)  are  likely  to  play  an  important  role  in  future  road 
traffic,  depending  on  the  predicted  future  regulation  scenarios  [1]. 
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Currently,  the  Lithium  ion  battery  technology  is  of  primary  interest 
to  the  automotive  industry,  since  it  provides  one  of  the  best  energy 
densities  available  today,  e.g.  ([2,3]). 

The  active  material  of  standard  Lithium  ion  battery  cells  is 
represented  by  stacked  or  wrapped  layers  of  the  cathode  and  anode 
material  sheets,  which  are  physically  separated  by  a  porous  but 
mechanically  robust  separator  foil,  in  order  to  help  prevent  an 
internal  short  circuit  of  the  cell.  Battery  systems  in  electric  vehicles 
are  protected  against  deformation  by  massive  structural  measures 
in  order  to  fulfil  federal  crash  laws,  consumer  crash  tests  and  other 
requirements.  Despite  these  achievements,  the  knowledge  of  the 
limiting  deformability  of  Lithium  ion  cells  prior  to  short  circuit 
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initiation  can  potentially  provide  valuable  additional  information 
for  the  design  of  future  battery  cells. 

On  the  experimental  side,  extensive  research  focused  on  the 
investigation  of  the  different  cell  components,  where  special 
emphasis  is  put  on  the  development  and  characterisation  of  sepa¬ 
rators  for  improved  cell  performance  ([4-6]). 

Simulation-based  research  included  stress  and  damage  analysis 
due  to  thermal  expansion  mismatch  of  the  different  electrode 
components  during  charging/discharging.  Xiao  et  al.  [7]  modelled 
the  stress  and  deformation  in  electrode  particles  and  separators 
using  a  meso-scale  representative  volume  element  (RVE)  coupled 
to  a  ID  macroscopic  battery  model.  Zhang  et  al.  [8]  simulated  the 
intercalation-induced  3D  stress  in  Lithium  ion  battery  electrode 
particles  of  different  size  and  shape.  They  concluded  that 
smaller  sized  particles  with  large  aspect  ratios  lead  to  reduced 
intercalation-induced  stresses  and  damage. 

Extensive  research  also  went  into  the  analysis  and  simulation  of 
abuse  tests  (oven,  short  circuit,  overcharge,  nail,  crush),  with 
emphasis  on  the  cell  behaviour  after  short  circuit  initiation.  A 
summary  paper  is  provided  by  Spotnitz  and  Franklin  [9].  Maleki 
and  Howard  [10]  investigated  the  internal  short  circuit  behaviour  of 
prismatic  cells,  subject  to  small  nail  penetration,  small  indentation 
and  cell  pinching,  experimentally  and  numerically.  They  concluded 
that  the  risk  of  a  critical  thermal  runaway  of  the  cell  is  mainly 
controlled  by  the  local  heat  conduction  of  the  cell  structure  around 
the  shorted  region.  In  the  worst  case  of  the  investigated  load  cases, 
namely  cell  pinching,  the  heat  cannot  be  transferred  to  the  cell-can, 
but  is  transferred  back  into  the  jelly  roll.  Spotnitz  et  al.  [11]  simu¬ 
lated  the  thermal  abuse  resistance  of  18  650-size  cells.  Analogous  to 
Maleki  and  Howard  [10],  the  authors  identified  heat  transfer  as  the 
governing  factor  for  thermal  runaway,  where  thermal  runaway  of 
a  single  cell  of  a  battery  pack  is  more  likely  to  cause  thermal 
runaway  of  the  whole  battery  pack  when  the  initiating  cell  is  in 
good  contact  with  other  cells  and  is  close  to  the  pack  wall. 

Safety  issues  related  to  Lithium  ion  battery  cells  are  discussed 
controversially  in  the  literature.  Generally,  there  are  two  different 
strategies. 

1.  Improving  the  safety  by  avoiding  abuse  and  protecting  the  cells 
appropriately. 

2.  Improving  the  safety  of  the  cells  by  understanding  the  mech¬ 
anisms  causing  internal  short  circuits  under  mechanical 
loading,  and  by  improving  the  cells  or  cell  components 
accordingly. 

Farrington  [12]  suggested  that  “effective  regulations  should 
promote  and  maximize  safe  transportation  of  lithium  batteries 
through  [...]  the  elimination  of  unsafe  circumstances”.  Sahraei 
et  al.  [13]  believe  that  “advanced  constitutive  models  are  needed 
for  strength/weight  optimization  and  safety  assessments  of  Li-ion 


batteries”.  The  authors  performed  a  comprehensive  structural 
testing  program  on  prismatic  pouch  cells  with  the  aim  to  set  up 
a  finite  element  model  for  representing  the  deformation  behaviour 
of  the  interacting  cell  components  on  a  macro-mechanical 
homogenized  scale,  and  also  on  a  refined  meso-mechanical  scale, 
where  individual  layers  of  the  electrodes  are  modelled  distinctively. 
In  another  paper  [14],  the  same  authors  investigated  cylindrical 
cells  subject  to  axial  crushing,  radial  crushing  (without  end  caps) 
and  radial  indentation,  and  used  numerical  finite  element  simula¬ 
tions  and  analytical  solutions  for  describing  the  deformation 
response  of  the  cells. 

The  two  papers  by  Sahraei  et  al.  ([13,14])  were  the  only  publi¬ 
cations  in  the  open  literature  found  by  the  present  authors  in  which 
the  deformation  behaviour  of  battery  cells  during  mechanical  abuse 
loading  was  investigated  and  modelled.  No  publications  could  be 
identified  for  predicting  the  onset  of  short  circuits  of  battery  cells 
under  mechanically  loading,  using  finite  element  modelling. 

In  the  present  study,  cylindrical  Lithium  ion  cells  are  subjected 
to  various  mechanical  abuse  tests,  in  order  to  create  a  crash 
simulation  model  of  the  cell.  The  developed  simulation  model 
allows  representation  of  the  cell  deformation,  and  also  features 
a  stress-based  fracture  criterion  for  predicting  the  load  state  and 
location  for  internal  short  circuit  onset  during  deformation. 

2.  Mechanical  cell  abuse  testing  program 

2.1.  Test  preparation 

Cylindrical  Nickel  Cobalt  Oxide  (NCA)  Lithium  ion  cells  (GAIA, 
HP  602030  NCA-45  Ah/162  Wh)  are  used  in  this  study.  The  main 
cell  dimensions  are  provided  in  Fig.  la,  and  the  interested  reader  is 
referred  to  the  online  data  sheet  of  the  cell  [15]  for  further  infor¬ 
mation  about  the  physical  and  mechanical  characteristics,  the 
electro-chemical  characteristics,  and  the  operating  conditions  of 
the  cell. 

The  cells  are  discharged  at  a  state  of  charge  (SoC)  of  slightly 
greater  than  0%,  or  to  the  recommended  voltage  limit  for  discharge, 
approximately  3  V,  respectively. 

A  zero  SoC  is  expected  to  enable  assessment  of  the  heat  location 
and  propagation  at  internal  short  circuit  initiation,  but  at  the  same 
time  allows  avoiding  severe  cell  reaction  like  smoke  and  fire 
development.  Three  temperature  sensors  are  applied  on  the  cell 
mandrel,  Fig.  la,  and  the  cell  voltage  is  also  recorded  during  the 
test. 

For  safety  reasons,  the  cells  are  located  in  a  windowed  cave 
featuring  a  gas  exhaust  device,  Fig.  lb. 

The  cave  is  mounted  on  a  hydraulic  tension/compression 
testing  device  allowing  the  application  of  compressive  loading 
using  different  types  of  punches  and  bearings,  Fig.  lc.  The  tests 
are  performed  under  quasi-static  loading  at  0.1  mm  s-1.  During 


Fig.  1.  a)  Main  cell  dimensions  (in  mm)  and  cell  instrumentation;  b)  test  chamber  set  up;  c)  investigated  load  cases. 
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Fig.  2.  Experimentally  obtained  punch  load-displacement  and  cell  voltage  curves  for 
the  load  case  “radial  crushing”  of  the  whole  cell. 
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Fig.  4.  Experimentally  obtained  punch  load-displacement  and  cell  voltage  curves  for 
the  load  case  “bending”  of  the  whole  cell. 


the  test,  the  load  and  machine  cross  head  movement  are  recor¬ 
ded.  An  infrared  camera  captures  the  front  view  of  the  cell  during 
the  test.  The  mechanical  test  was  interrupted,  once  the  cell 
voltage  dropped  to  zero,  indicating  a  potential  internal  short 
circuit.  At  this  state,  the  punch  was  removed  from  the  cell,  and 
the  temperature  development  of  the  cell  was  monitored  another 
hour  after  the  test,  in  order  to  prove  the  occurrence  of  an  internal 
short  circuit.  Each  load  case  configuration  has  been  repeated  at 
least  two  times. 


2.2.  Test  analysis 

The  load-displacement  curves  and  the  cell  voltage 
measurements  are  provided  in  Figs.  2-4.  It  is  noted  that  the 
cells  could  endure  significant  deformation  prior  to  the  occurrence 
of  an  internal  short  circuit.  A  good  repeatability  of  the 
load-displacement  curves  and  of  the  punch  displacement,  at 
which  the  cell  voltage  drops  significantly,  is  observed.  In  some 
cases  the  cell  voltage  did  not  immediately  decrease  to  zero,  but 
remained  at  a  reduced  level.  In  these  cases,  the  unloaded  punch  has 
been  applied  again  and  driven  further  until  the  voltage  signal 
decreased  to  zero,  Fig.  2:  Expl.  In  one  of  the  bending  tests,  unde¬ 
sired  movement  of  the  clamped  bearings  has  been  observed,  Fig.  4: 
Exp3. 

For  all  tests,  the  temperature  sensors  recorded  an  increasing 
temperature  from  about  25  °C  to  approximately  50  °C  after  the 
tests,  proving  the  occurrence  of  an  internal  short  circuit.  Addi¬ 
tionally,  the  infrared  camera  recorded  high  temperature  spots  at 
the  beginning  of  the  short  circuit  (cell  voltage  drop),  Fig.  14  in 
Section  3.5,  revealing  the  location  of  the  short  circuit.  This  crucial 
information  is  used  as  a  quality  check  of  the  numerical  short  circuit 
criterion  developed  in  Section  3.5. 


Fig.  3.  Experimentally  obtained  punch  load-displacement  and  cell  voltage  curves  for 
the  load  case  “indentation”  of  the  whole  cell. 


Interestingly,  a  sudden  drop  in  the  force  level  is  observed  in 
some  of  the  load-displacement  curves,  just  before  initiation  of 
the  internal  short  circuit.  The  effect  is  most  pronounced  for  the 
load  case  “indentation”,  Fig.  3.  The  steel  housing  did  not  show 
any  sign  of  fracture.  Hence,  it  is  assumed  that  the  load  drop  is 
caused  by  structural  fracture  within  the  jelly  roll.  A  computer 
tomography  (CT)  analysis  of  the  cell  subjected  to  bending  loading 
indeed  revealed  a  macroscopic  fracture  of  the  jelly  roll,  Fig.  5.  It  is 
interesting  to  recapitulate  that  the  onset  of  the  short  circuit  is 
not  just  due  to  fracture  phenomena  of  the  separator  on  a  micro- 
or  meso-scale.  Macroscopic  through-thickness  fracture  of  the 
jelly  roll  enables  relative  movement  of  the  “unprotected”  anode 
and  cathode  material  at  the  fracture  surface,  which  then  easily 
can  lead  to  direct  physical  contact  of  the  anode  and  cathode,  or 
to  bridging  contact  via  the  steel  housing,  initiating  the  short 
circuit. 

For  the  other  load  cases,  where  the  zone  of  internal  short 
circuit  is  subject  to  compressive/shear  loading,  no  such  clear 
statement  can  be  made  regarding  the  existence  of  a  macroscopic 
fracture.  It  is  however  likely  that  existing  cracks  become  closed 
due  to  the  compressive  state  of  loading,  and  hence  are  difficult  to 
detect. 

In  the  next  section,  a  numerical  short  circuit  detection  criterion 
is  developed,  based  on  the  assumption  that  fracture  within  the  jelly 
roll  initiates  the  short  circuit. 


3.  Deformation  and  fracture  modelling  of  the  cell 

3.1.  Preparation  of  the  simulation  model 

The  commercial  explicit  crash  code  Virtual  Performance  Solu¬ 
tion  (VPS)  is  used  for  the  numerical  analysis.  An  overview  of  the 
finite  element  model  of  the  cell,  punches  and  bearings  is  provided 
in  Fig.  6. 

The  jelly  roll  is  secured  by  a  1  mm  thick  cylindrical  stainless 
steel  housing.  The  housing  is  represented  by  shell  elements 


Fig.  5.  Deformed  cell  and  computer  tomography  for  the  load  case  “bending”. 


380 


L.  Greve,  C.  Fehrenbach  /  Journal  of  Power  Sources  214  (2012)  377-385 


Component 

Thickness 

Material 

Plasticity  model 

Remarks 

Cap 

0.2  mm 

Stainless  steel 

Isotropic  J2  yield  criterion 

No  influence  on 
simulation 
results;  just 
included  for 
optical  reasons. 

Outer  skin  1.0  mm  Stainless  steel  Hill48  yield  criterion  with  Shell  elements 

normal  anisotropy  with  5  through¬ 

thickness 
integration 
points 


Jelly  roll 

Nickel  Cobalt 

Pressure-dependent 

Central  hole  and 

Oxide  (NCO) 

isotropic  Raghava  yield 

pole 

electro- 

criterion  with  non- 

connections  are 

chemestry 

associated  flow  rule 

neglected. 

Fig.  6.  Macro-mechanical  finite  element  model  of  the  cell  (using  quarter  model  due  to  symmetry  conditions). 


(approximately  5x5  mm2),  whereas  the  jelly  roll  is  represented  by 
volume  elements  (approximately  5x5x5  mm3).  The  inner  hole  of 
the  jelly  roll  is  neglected  in  the  simulation  model,  as  well  as  dis¬ 
cretization  of  the  pole  connectors  on  each  end  of  the  jelly  roll.  The 
bearings  and  the  punches  are  modelled  as  rigid  bodies.  A  penalty¬ 
spring  based  contact  definition  is  applied  between  housing  and 
jelly  roll,  assuming  a  Coulomb  friction  coefficient  of  0.1.  Another 
contact  is  applied  between  the  punch,  the  cell  housing  and  the 
bearings,  using  a  friction  coefficient  of  0.1. 

The  explicit  simulations  were  performed  using  a  time  step  of 
one  micro-second,  which  corresponds  to  a  typical  value  used  for 
full  scale  vehicle  crash  simulations.  The  density  of  the  jelly  roll  was 
then  trimmed  in  order  to  represent  the  real  cell  weight  of  1.5  kg  in 
the  simulation  run.  The  punch  intrusion  velocity  was  set  to 
a  significantly  higher  value  (vsim  =  0.5  mm  ms-1)  than  the  actual 
test  speed  (v£Xp  =  0.1  mm  s_1),  in  order  to  reduce  the  computation 
time.  At  any  time  during  the  simulation,  the  difference  of  punch 
force  and  bearings  force  was  less  than  0.2  kN  (0.04%  of  the 
maximum  load),  indicating  that  existing  mass  inertia  forces  in  the 
simulation  are  still  negligible  in  the  simulated  configuration.  All 
material  models  used  in  the  subsequent  sections  are  rate  inde¬ 
pendent,  so  that  the  increased  simulation  speed  does  not  influence 
the  model  behaviour  at  all. 

3.2.  Identification  of  material  model  parameters  for  the  housing 

Within  Volkswagen  Group  Research,  a  modular  user  material 
model  (MMM)  for  the  crash  simulation  of  deformation  and  fracture 
of  automotive  materials  has  been  developed  in  the  past.  The 


material  model  is  based  on  a  modular  concept,  where  the  user  can 
plug-in  different  modules  for  elasticity,  yield  criterion,  hardening, 
fracture  and  damaging,  as  required  for  the  given  material  of 
investigation.  In  this  study,  MMM  is  used  for  the  simulation  of  all 
cell  components. 

Tensile  specimens  are  cut  from  flattened  housing  sheets,  Fig.  7. 
The  tensile  tests  were  performed  on  an  electromechanical  tensile 
machine  at  0.01  mm  s-1,  where  the  load  and  the  local  displacement 
within  the  gauge  section  were  recorded  during  the  test.  Textbook 
values  are  used  for  the  Young’s  modulus,  Ehousing  =  207  GPa  and  the 
elastic  Poisson  ratio,  v  =  0.3. 

The  tensile  tests  revealed  a  normal  anisotropy  of  the  housing 
sheet,  r  =  0.64,  which  is  represented  by  the  classical  Hill48  yield 
criterion  [16].  The  strain  hardening  could  be  well  represented  by 
the  modified  Swift  law  [17],  as  a  function  of  the  equivalent  plastic 
strain  W, 

housing (^)  -  (1) 

with  k  =  0.7  GPa,  r0  =  0.00801  and  the  hardening  exponent 
n  =  0.1385,  Fig.  7b. 

A  three-point  tube  bending  test  is  performed  for  validation  of 
the  identified  material  model  parameters  of  the  housing,  showing 
very  good  agreement  of  the  simulation  and  experiment,  Fig.  8. 

3.3.  Constitutive  model  for  the  jelly  roll 

The  jelly  roll  is  represented  by  a  wrapped  stacked  sequence  of 
separator,  anode,  separator  and  cathode  foils,  Fig.  9a  and  b.  The 


Fig.  7.  a)  Specimen  extraction  from  the  flattened  housing  tube;  b)  experimental  and  analytical  hardening  curve. 
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Fig.  8.  Three-point  bending  test  of  cylindrical  housing  tube:  comparison  of  experi¬ 
mental  and  simulated  load-displacement  curves. 


anode  material  (LiCo02  powder)  is  baked  onto  the  surfaces  of  the 
Aluminium  current  collector,  whereas  the  cathode  material 
features  a  central  Copper  foil  current  collector  with  a  graphite 
powder  coating,  Fig.  9b.  Both  electrode  coating  materials  are 
highly  porous,  and  soaked  in  electrolyte  fluid.  A  detailed 
numerical  analysis  would  require  a  comprehensive  character¬ 
isation  of  all  jelly  roll  components  (cathode,  anode,  separator  and 
foils),  and  their  interaction,  e.g.  friction  and  inter-laminar  slip. 
Since  the  thickness  of  a  separator  foil  is  typically  about  0.03  mm, 
a  full  3D  finite  element  representation  of  the  unit  cell  would 
require  finite  elements  featuring  an  element  length  of  similar 
geometrical  size,  eventually  leading  to  an  unmanageable  number 
of  finite  elements  and  computational  requirements  for  industrial 
crash  simulation  application.  Hence,  in  this  study  a  simplified 
macro-mechanical  approach  is  applied,  in  which  the  jelly  roll  is 
considered  isotropic,  Fig.  9c.  That  is,  the  deformation  of  the 
electrodes  together  with  the  separator  is  treated  in  an  integrated 
and  continuous  way.  The  influence  of  frictional  slip  between  the 
foils  and  the  influence  of  potential  foil  separation  (gaping)  is 
neglected. 

During  compressive  mechanical  testing,  the  electrolyte  is 
squeezed  out  of  the  porous  electrodes,  and  moves  into  the 
empty  cavitation  at  the  ends  of  the  cells,  so  the  porous  electrodes 
become  compressed  and  change  their  circumscribing  volume.  For 
phenomenological  modelling  of  the  electrode  compressibility, 
a  pressure  dependent  yield  criterion,  and  a  pressure  dependent 
non-associated  flow  rule  are  proposed.  In  the  following  procedure, 
a  generalized  plasticity  model  for  the  jelly  roll  components  will  be 
proposed.  Then,  the  generalized  model  will  be  simplified  for  the 
present  isotropic  modelling  approach  according  to  Fig.  9c.  The  yield 
criterion  reads: 


Gray-shaded  zone: 
Jelly  roll 


a  b 


Fig.  10.  Analysis  of  the  axial  jelly  roll  deformation  at  the  initial  and  final  state:  a)  CT 
analysis:  b)  jelly  roll  deformation  obtained  by  simulation. 


/  =  a  -  ^jellyroll(cp)  =  0.  (2) 

where  ftjeiiyroii  is  the  hardening  function  of  the  jelly  roll,  and  a  is 
the  equivalent  stress,  represented  by  a  modified  Raghava  yield 
criterion  [18]: 


a  =  2^(3(“- 1)(Tm  +  \/9(“-1)2<+4«  <^ev),  (3) 

in  which  the  original  von  Mises  equivalent  stress  is  replaced  by 
a  generalized  anisotropic  deviatoric  equivalent  stress  adev. 

In  equation  (3),  a  represents  the  absolute  ratio  of  uniaxial 
compressive  and  uniaxial  tensile  yield  strength  of  the  material,  and 
am  is  the  hydrostatic  stress, 


Om  =  3^1  +  v2  +  <73), 


(4) 


where  Oi  are  the  principal  stresses.  adev  represents  an  anisotropic 
quadratic  yield  criterion,  obtained  by  linear  transformation  of  the 
stress  vector, 


^dev  =  l[P(a-b)]r(d-b),  (5) 

where  P  is  the  linear  transformation  matrix, 
a  =  [au  g2 2  033  <r12  023  0"i3  ]T  is  the  stress  vector,  and  b  = 
[bu  b22  b33  0  0  0]T  is  the  initial  back  stress  vector. 


Fig.  9.  a)  Cross  section  of  the  jelly  roll;  b)  unit  cell  of  the  jelly  roll  components;  c)  simplified  isotropic  continuum  representation  of  the  jelly  roll  material. 
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With  an  appropriate  choice  of  the  matrix  components  of  P,  the 
deviatoric  equivalent  stress  represents  the  classical  Hill48  yield 
criterion  with  an  initial  back  stress  vector, 


corresponds  to  a  low  plastic  Poisson’s  ratio  in  the  material  model.  A 
non-zero  value  of  *^=0.15  is  used  for  numerical  stability  in  the 
simulation  model. 


^dev  — 


\j^{a22b  ~  a33b)2+G(a33b  ~  Gnb)2+H(aUb  ~  G22 b)2+2  [^(^23 )2 (^1 3 )2 +^(^1 2 )2]  5 


(6) 


with 


°\\b  =  °\\  ~  ^li;  a22b  =  G22~b22',  a33b  =  °33  ~  ^33,  C7) 

The  Hill  parameters  (F,  G,  H,  L,  M,  N )  would  have  to  be  identified 
from  experiments. 

The  proposed  anisotropic  model  provides  a  flexible  general 
model  for  representation  of  the  yield  behaviour  of  the  electrode 
particles  and  separator  materials  on  a  meso-modelling  scale. 
However,  identification  of  all  material  parameters  would  require 
extensive  testing,  which  is  not  scope  of  the  present  macro¬ 
mechanical  study.  For  the  subsequent  analysis,  the  general  model 
is  significantly  simplified.  Due  to  the  assumption  of  isotropy,  the 
Hill  parameters  are  set  to  their  isotropic  values,  F=G  =  H=  0.5,  and 
L  =  M  =  N  =  1.5,  and  the  back  stress  vector  is  b  =  0.  No  data  was 
available  for  the  determination  of  a  potential  yield  strength 
difference,  hence  a  =  1.  So  eventually,  the  yield  criterion,  equation 
(3),  for  the  present  macro-mechanical  isotropic  modelling 
approach  is  simplified  to  standard  J2  (von  Mises)  theory, 


For  the  investigated  load  cases,  the  jelly  roll  deformation  occurs 
mainly  in  radial  compression  and  in  radial  shear  direction. 
Considering  the  underlying  monotonic  mechanical  loading,  the 
differentiation  of  elastic  and  plastic  strains  is  not  necessary,  but  for 
the  elasto-plastic  constitutive  model  featuring  additive  strain 
decomposition,  the  elastic  properties  have  to  be  input  separately. 
Due  to  the  relatively  large  Young’s  modulus  of  the  current  collector 
foils,  the  jelly  roll  deformation  is  mainly  taking  place  in  the  non- 
metallic  porous  particle  coatings  of  the  jelly  roll.  Hence,  a  small 
Young’s  modulus  is  chosen  for  the  jelly  roll,  using 
Ejeiiyroii  =  1.5  GPa,  which  is  calibrated  along  with  the  plastic  strain 
hardening  curve,  for  which  the  following  analytical  function  is 
proposed: 


R(lP)  = 


^plateau  (^plateau  ^yield)  x  Exp 
X  [l+s(^)m]. 


^ref 


(12) 


-ff22)2  +  (ffll  -ff33)2  +  (^22 -ff33)2+6[(ff23)2  +  (<713)2  +  (<712)2]}- 


(8) 


A  non-associated  pressure  dependent  flow  rule  is  proposed, 

As p  =  AA— ,  (9) 

CT 

where  Ad3  is  the  plastic  strain  increment  vector,  AA  is  the  plastic 
multiplier,  and  g  is  the  flow  rule  function, 

S  =  \/ffdev  +  02ffm-  (10) 

Shape  parameter  can  be  related  to  the  plastic  Poisson  ratio  vp 
under  uniaxial  tension, 


The  first  term  represents  the  first  phase  of  the  hardening 
behaviour,  from  initial  yield  <7yieid,  up  to  the  intermediate  plateau 
stress,  0-piateau,  where  rref  shapes  the  curvature  of  the  hardening 
curve  in  between.  The  second  term  represents  subsequent 
progressive  strain  hardening,  controlled  by  the  power  law  param¬ 
eters  s  and  m.  The  calibrated  parameters  are  cryieid  =  0.0003  GPa, 
o plateau  =  0.0006  GPa,  £ref  =  0.009,  s  =  1000  and  m  =  2.7. 

For  the  load  case  “radial  crushing”  the  computed  cell  deforma¬ 
tion  is  provided  in  Fig.  10b,  being  in  good  agreement  with 
the  experiment.  The  comparison  of  the  load-displacement 
curves  is  provided  in  Fig.  11a,  accompanied  with  the  predicted 
load-displacement  curves  for  the  load  cases  “bending”  (Fig.  lib) 
and  “indentation”  (Fig.  11c).  All  simulation  curves  are  in  good 
agreement  with  the  corresponding  experiment. 


P 


3_ 

71 


1-2  vP 

TT ~W' 


(ID 


3.5.  Development  of  a  fracture  and  short  circuit  criterion  for  the 
jelly  roll 


3.4.  Parameter  identification  for  the  jelly  roll,  and  simulation  of  the 
cell  tests 

The  radial  crushing  load  case  is  used  for  calibration  of  the 
plasticity  model.  CT  analysis  revealed  that  whilst  the  jelly  roll 
becomes  significantly  compacted  during  the  test,  it  exhibited 
almost  no  axial  straining.  The  distance  from  the  jelly  roll  to  the 
boundary  of  the  housing  remained  roughly  the  same,  Fig.  10a.  That 
is,  the  jelly  roll  exhibits  almost  perfect  compressibility,  which 


As  discussed  in  Section  2.2,  macro-mechanical  fracture  of  the 
jelly  roll  is  clearly  revealed  from  the  CT  analysis  for  the  load  case 
“bending”,  and  also  indicated  by  a  significant  sudden  decrease  of 
the  punch  force  followed  by  the  initiation  of  an  internal  short 
circuit  for  the  load  case  “indentation”,  Fig.  3.  It  is  assumed  that 
fracture  of  the  jelly  roll  leads  to  internal  short  circuit  in  all  inves¬ 
tigated  load  cases,  since  once  fracture  is  initiated,  the  relative 
movement  of  the  separated  jelly  roll  parts  easily  leads  to  a  short 
circuit  initiating  connection  of  anode  and  cathode  material. 
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Fig.  11.  Comparison  of  experimental  and  computed  load-displacement  curves,  and 
the  recorded  voltage-displacement  for  the  different  load  cases:  a)  “radial  crushing”: 
b)  “bending”:  c)  “indentation”. 


In  the  fracture  zone,  the  jelly  roll  is  subject  to  compression- 
shear  loading  for  the  load  cases  “radial  crushing”  and  “indenta¬ 
tion”  and  subject  to  tensile  loading  for  the  load  case  “bending”. 
Hence,  a  fracture  criterion  is  sought,  which  can  take  into  account 
potential  differences  in  the  fracture  behaviour  under  compressive/ 
shear  and  normal  stresses.  The  classical  stress-based  fracture 
criterion  by  Coulomb  [19]  and  Mohr  [20]  can  describe  such 


Fig.  12.  Mohr-Coulomb  fracture  criterion:  a)  shear  and  normal  stress  acting  on  an 
arbitrary  plane  in  principal  stress  space;  b)  graphical  visualization  of  the  solved 
maximum  value  problem  in  principal  stress  space,  and  the  definition  of  scaling 
distance  A. 


radial 

bend 

indent 

■MC 


cr3  [GPa] 

Fig.  13.  Calibration  procedure  for  the  MC  fracture/short  circuit  criterion. 


behaviour.  Hence,  the  criterion  is  reviewed  and  applied  subse¬ 
quently.  The  Mohr-Coulomb  (MC)  criterion  for  an  isotropic 
continuum  reads 

maxJx  +  CTff,,}  =  c2,  (13) 

where  t  and  an  are  the  shear  and  normal  stresses  acting  on  a  cutting 
plane  in  principal  stress  space,  Fig.  12a.  Fracture  occurs  on  the  plane 
perpendicular  to  the  normal  vector  ~n,  where  the  shear  stress  t  plus 
the  normal  stress  on  scaled  by  C\  represents  the  maximum  C2.  In 
equation  (13),  C\  and  c2  are  material  constants,  which  have  to  be 
calibrated  from  experiments. 

The  maximum  value  problem  can  be  solved  analytically,  which 
has  been  done  by  Bai  and  Wierzbicki  [21],  for  instance,  leading  to 
the  simple  expression 


°T  =  P&3  +  2c2,  (14) 

where  o\  and  a3  are  the  first  and  third  principal  stresses,  respec¬ 
tively,  and 


p  =  (\d+ci  ~ci) 
(\A +ci +ci) 


(15) 


Equation  (15)  represents  the  fracture  line,  where  fracture  occurs 
once  the  combination  of  first  and  third  principal  stresses  exceeds 
the  line,  Fig.  12b.  For  the  simulation  model,  a  scaling  distance  A  is 
introduced, 


(16) 


which  represents  the  scaling  distance  from  the  current  stress  state 
~a,  to  the  fracture  (short  circuit)  line,  ~o\  assuming  a  constant  stress 
ratio,  Fig.  12b.  Jelly  roll  fracture,  accompanied  by  an  internal  short 
circuit,  occurs  at  A  >  1. 

The  calibration  procedure  for  the  MC  model  parameters  C\  and 
c2  is  as  follows: 

•  Running  the  simulation  models  for  all  load  cases  to  the  punch 
displacement  at  which  a  short  circuit  is  observed  in  the 
experiment. 

•  Visualize  the  final  principal  stresses  (o\  and  o3)  of  all  finite 
elements  of  the  jelly  roll  in  one  chart,  Fig.  13. 

•  Determine  the  fracture  line  by  calibrating  the  MC  parameters 
(Cl,  c2),  Fig.  13. 


The  obtained  values  are  Ci  =  0.865  and  c2  =  0.012  GPa. 
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Load  case  Short  (Exp)  Short  (Sim) 


Fig.  14.  Comparison  of  the  experimentally  observed  short  circuit  locations  and  corresponding  model  predictions  for  the  different  load  cases:  a)  “radial  crushing”:  b)  “indentation”: 
c)  “bending”. 


It  is  noted  that  the  calibration  of  the  MC  criterion  is  solely  based 
on  stress  analysis,  where  the  location  of  the  short  circuit  has  not 
been  taken  into  account  by  any  means.  Hence,  it  is  now  interesting 
to  compare  the  predicted  short  circuit  locations  with  the  experi¬ 
mentally  observed  IR  measurements  provided  in  Fig.  14,  in  which 
high  and  highest  temperatures  are  coloured  red  and  white, 
respectively. 

The  experimental  infrared  picture  of  the  load  case  “radial 
crushing”,  Fig.  14a,  indicates  short  circuit  initiation  at  one  end  of  the 
jelly  roll,  close  to  the  pole  connection.  Examination  of  the  longi¬ 
tudinal  sectional  plane  of  the  corresponding  simulation  shown  in 
Fig.  14a  reveals  that  in  the  core  zone  of  the  jelly  roll,  the  state  of  A  is 
quite  homogeneous  along  the  cell  axis,  and  close  to  the  limiting 
value  of  unity,  at  which  fracture  and  internal  short  circuit  is  initi¬ 
ated.  Therefore,  it  is  assumed  that  the  location  of  the  short  circuit 
initiation  is  eventually  triggered  by  some  mechanical  interaction  of 
the  pole  connector  ant  the  end  of  the  jelly  roll. 

For  the  load  case  “indentation”  the  short  circuit  location  is 
predicted  close  to  the  punch,  whereas  for  the  load  case  “bending” 
the  location  is  predicted  on  the  surface  of  the  jelly  roll  opposite  to 
the  punch,  Fig.  14b  and  c.  Both  predictions  are  in  good  agreement 
with  the  experimental  observations,  which  indicates  that  the 
applied  MC  criterion  provides  a  simple  but  useful  method  for  the 
detection  of  jelly  roll  fracture  and  subsequent  short  circuit 
initiation. 


4.  Conclusions 

A  quasi-static  mechanical  test  program  on  commercial  cylin¬ 
drical  Lithium  ion  cells  has  been  carried  out,  in  order  to  assess  the 
deformability  until  initiation  of  an  internal  short  circuit.  A  macro¬ 
mechanical  finite  element  model  of  the  cell,  suitable  for  standard 
crash  simulations,  has  been  set  up,  using  an  element  mesh  size  of 
approximately  5  mm.  The  material  model  properties  for  the  steel 
housing  have  been  obtained  from  tensile  tests,  and  the  material 
properties  for  the  jelly  roll  have  been  calibrated  from  one  of  the  cell 
tests,  namely  the  radial  crushing  load  case.  The  classical  stress- 
based  fracture  criterion  after  Mohr  and  Coulomb  has  been  uti¬ 
lised  for  fracture  modelling  of  the  jelly  roll,  which  is  considered 
isotropic  in  the  present  study.  It  is  shown  that  the  Mohr— Coulomb 
criterion  correctly  predicts  both  the  punch  displacement  to  frac¬ 
ture,  and  also  the  correct  different  fracture  locations  for  every  load 


case.  The  predicted  fracture  locations  correspond  to  the  observed 
locations  of  the  internal  short  circuits  of  the  cells.  The  testing 
demonstrated  that  the  cells  could  withstand  considerable  defor¬ 
mation  prior  to  the  onset  of  an  internal  short  circuit. 

Future  research  should  involve  the  experimental  character¬ 
isation  and  numerical  simulation  of  the  individual  jelly  roll 
components  (anode,  separator  and  cathode)  and  their  interaction 
during  mechanical  loading.  Fracture  analysis  of  the  separator  on 
a  meso-scale  could  inspire  the  development  of  improved  macro¬ 
scale  models  and  criteria  suitable  for  industrial  applications. 
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